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Precipitate hardening is commonly used in materials science to control strength by acting on the number den- 
sity, size distribution, and shape of solute precipitates in the hardened matrix. The Fe-Cu system has attracted 
much attention over the last several decades due to its technological importance as a model alloy for Cu steels. 
In spite of these efforts several aspects of its phase diagram remain unexplained. Here we use atomistic simu- 
lations to characterize the polymorphic phase diagram of Cu precipitates in body-centered cubic (BCC) Fe and 
establish a consistent link between their thermodynamic and mechanical properties in terms of thermal stability, 
shape, and strength. The size at which Cu precipitates transform from BCC to a close-packed 9R structure is 
found to be strongly temperature dependent, ranging from approximately 4 nm in diameter (^ 2, 700 atoms) 
at 200 K to about 8 nm (^ 22, 800 atoms) at 700 K. These numbers are in very good agreement with the inter- 
pretation of experimental data given Monzen et al [Phil. Mag. A 80, 711 (2000)]. The strong temperature 
dependence originates from the entropic stabilization of BCC Cu, which is mechanically unstable as a bulk 
phase. While at high temperatures the transition exhibits first-order characteristics, the hysteresis, and thus the 
nucleation barrier, vanish at temperatures below approximately 300 K. This behavior is explained in terms of 
the mutual cancellation of the energy differences between core and shell (wetting layer) regions of BCC and 9R 
nanoprecipitates, respectively. The proposed mechanism is not specific for the Fe-Cu system but could gener- 
ally be observed in immiscible systems, whenever the minority component is unstable in the lattice structure of 
the host matrix. Finally, we also study the interaction of precipitates with screw dislocations as a function of 
both structure and orientation. The results provide a coherent picture of precipitate strength that unifies previous 
calculations and experimental observations. 

PACS numbers: 62.20.fq, 64.70.kd, 64.70.Nd, 81.30.Hd 



I. INTRODUCTION 

Copper impurities contribute to high-temperature embrittle- 
ment in reactor pressure vessels (RPV) ferritic steels, which is 
a serious concern regarding lifetime extension of existing nu- 
clear reactors.^ This behavior originates from the very small 
solubility of Cu in a-Fe (body-centered cubic, BCC) and the 
correspondingly large driving force for precipitation. Cop- 
per has been observed to initially form small vacancy-rich 
clusters, which dissolve over time under typical RPV condi- 
tions, leaving embryonic Cu clusters that interact with internal 
dislocations.^ The evolution of these nano- sized precipitates 
must be taken into account in embrittlement models for life- 
time predictions of RPV steels.!^ 

Binary Fe-Cu alloys can be used to understand Cu precip- 
itation and dislocation-precipitate interactions as prototypes 
for more complex RPV steels.^ Based on a series of careful ex- 
perimental studies on these alloys, three transformation stages 
have been identified^ ^ : (/) Cu precipitates nucleate in the BCC 
structure, (//) as they grow they undergo a martensitic phase 
transition to a multiply-twinned 9R structure, which (///) even- 
tually transforms — presumably via diffusion — into the 3R 
structure.^ This sequence was determined using post-mortem 
transmission electron microscopy (TEM) investigations and 
selected area diffraction .'^'^'^^ 

Several molecular dynamics (MD) and statics simulation 
studies ^^^^ have been carried out to investigate the structural 
transitions of precipitates and to shed light on the underly- 
ing atomic mechanisms as such information is difficult to 



access via post-mortem experimental probes. A number of 
these studies have documented the shear instability associated 
with the aforementioned phase transform ation that develops 
in BCC precipitates above a certain size.G^™^ j^ j^^^ ^j^q 
been shown that vacancies bind strongly to copper precipi- 
tates, affect their stabilit>^^ ^^ and play a role in the loss of co- 
herency between precipitate and matrix^E^C^. Furthermore, 
motivated by the importance of Cu precipitates for the me- 
chanical properties of a-iron several investigations have ad- 
dressed their interaction with dislocations, in particular in the 
case of coherent BCC precipitates.^"^ ^^ ^^ 

The objective of the present work is to obtain a comprehen- 
sive picture of the correlation between structure and size of Cu 
precipitates in a-iron as a function of temperature, and to elu- 
cidate the interaction of dislocations with 9R precipitates. To 
this end, we first obtain the size-temperature phase diagram 
of Cu precipitates using a massively parallel hybrid molecular 
dynamics/Monte Carlo (MD/MC) algorithm.^^ This investiga- 
tion reveals a strong temperature dependence of the precipitate 
size at which the BCC-9R transition occurs, which originates 
from the entropic stabilization of BCC Cu. At high temper- 
atures the character of the transformation is decidedly first- 
order as it exhibits a pronounced hysteresis and a large latent 
heat. In contrast, at low temperatures both hysteresis and nu- 
cleation barrier vanish. At a microscopic scale the transition, 
however, appears to remain discontinuous. 

Regarding the dislocation-precipitate interaction, we find 
that the critical shear strength required for a dislocation to 
pass a precipitate changes discontinuously across the BCC- 



9R transition by approximately a factor of two and exhibits a 
marked dependence on orientation on the 9R side of the tran- 
sition. These results have important implications for the inter- 
pretation of post-mortem experiments carried out near room 
temperature in connection with material structure and perfor- 
mance under RPV operating conditions at temperatures in ex- 
cess of 560 K.^ 

The paper is organized as follows. Section |Il| describes our 
simulation methodology. Precipitate structure and transition 



mechanism are described in Sect. Ill The thermodynamics of 
the transition are detailed in Sect. IV followed by an analy- 
sis of the interaction of screw dislocations with BCC and 9R 
precipitates in Sect. |V| Finally, in Sect. VI the results are dis- 
cussed with respect to thermodynamic implications, their im- 
pact on the interpretation of experimental data, and plasticity. 



II. METHODOLOGY 

The Fe-Cu binary alloy system was modeled using 
semi-empirical interatomic potentials of the embedded atom 
method type. For Cu and Fe potentials by Mishin et al^ 
and Mendelev et al (version v2 in Ref. 23 ) were used, re- 
spectively, while Fe-Cu interactions were described using the 
potential by Pasianot and Malerba.^"^ This combination of po- 
tentials provides an accurate description of the alloy phase di- 
agram and yields very good agreement with density functional 
theory calculations for the energetics of small Cu clusters in 
BCC Fe.ISl 

Simulations were carried out using the hybrid MD/MC al- 
gorithm introduced in Refs. 21 as implemented in the mas- 
sively parallel MD code lammps.^^ The approach is based 
on the variance-constrained semi-grand-canonical (VC-SGC) 
ensemble,^ ^ ^^ which generalizes the semi-grand-canonical 
(SGC) ensemble to allow stabilization of multiphase equilib- 
ria. Orthorhombic simulation cells were employed with cell 
vectors oriented along [111], [110], and [112]. The majority 
of the simulations was based on simulation cells comprising 
42 X 52 X 45 repetitions of the unit cells corresponding to 
786,240 atoms with additional simulations of smaller precip- 
itate sizes at low temperatures using 31 x 38 x 33 unit cells 
(310,992 atoms). Periodic boundary conditions were applied 
in all directions. The Nose-Hoover thermostat and barostat^^ 
were employed to equilibrate the system at ambient pressure 
and temperatures between 300 and 700 K. The equations of 
motion were integrated using a timestep of 2.5 fs. Simulations 
were run for at least 1.6 x 10^ MD steps and — if necessary 
for reaching equilibrium and/or collecting more statistics — 
were extended to up to 4 x 10^ MD steps. Every 20 steps 
the MD simulation was interrupted to carry out a number of 
VC-SGC MC trial moves corresponding to 10% of a full MC 
sweep. Each simulation run thus comprises between 8,000 
and 20,000 attempts per particle to swap the atom type. 

Parameters for the VC-SGC MC simulations were obtained 
as follows. First a series of SGC-MC simulations was carried 
out to determine the relation between relative chemical poten- 
tial A/i and concentration at temperatures between 300 and 
700 K. The average constraint parameter (p in the VC-SGC 



MC scheme^ ^ was then taken as the A/i value at which the 
solid solution of Cu in Fe became unstable. The variance con- 
straint parameter k, was set to 2000 eV/atom/ZcB^.^^ Finally, 
the target concentration cq was varied to obtain the desired to- 
tal number of Cu atoms in the simulation cell and thereby the 
size of the Cu precipitate. 

Simulations of dislocation-precipitate interactions were 
conducted at a temperature of 300 K with open boundary con- 
ditions perpendicular to the shear plane. Shear stress was ap- 
plied so as to produce a Peach-Kohler force on the (IIO)bcc 
plane that intersects the Cu precipitate through its equatorial 
plane. 

Structural changes in the simulations were detected using 
the Ackland- Jones parameter,^^ which uses bond angle distri- 
butions to classify local environments as body-centered cubic 
(BCC), face-centered cubic (FCC), hexagonal-close-packed 
(HCP), or icosahedral. The reliability of the Ackland- Jones 
analysis was enhanced by position- averaging over blocks of 
800 MD steps. Precipitate structures were analyzed and visu- 
alized using OVITOP 



III. STRUCTURAL ASPECTS OF BCC-9R TRANSITION 
A. General features of BCC and 9R precipitates 

Two different types of precipitates were observed in the 
simulations, examples for which are shown in Fig.[T] The pre- 
cipitate shown in Fig. [TJb) represents a snapshot from a VC- 
SGC MC/MD simulation at 700 K, in which the total concen- 
tration was fixed at 2%, a value significantly above the bulk 
solubility.^ ^ This overs aturation leads to formation of a spher- 
ical Cu precipitate containing approximately 14,900 atoms, 
which is practically Fe free. On average the precipitate is 
found to retain the BCC structure of the surrounding matrix. 

The precipitate shown in Fig. [TJe) was also obtained at 
TOOK using, however, a target concentration of 4.5%. It ex- 
hibits an elongated ellipsoidal shape, contains about 34,600 
atoms and appears to be composed of ordered stacks of FCC 
and HCP atoms, which can be identified as belonging to a 
multiply-twinned 9R structure as illustrated in Fig. [l|g)p2 
There are three distinct ways in which close-packed planes 
can be stacked denoted here by A, B, and C. For exam- 
ple, the FCC structure is described by the stacking sequence 
ABC — ABC whereas the sequence AB — AB yields the HCP 
structure. By contrast the unit cell of the 9R structure com- 
prises nine close-packed planes. Since the Ackland- Jones pa- 
rameter considers only nearest neighbors it is suitable for dis- 
tinguishing HCP and FCC environments but cannot directly 
identify the 9R structure. As shown in Fig. [TJg) the 9R struc- 
ture rather appears as a sequence of FCC and HCP atoms. 
Comparison of the ideal 9R structure in Fig. [TJg) with the 
precipitates in Fig. [TJd-f) demonstrates that the latter have 
adopted a multiply-twinned 9R structure. 

It is important to note that 9R precipitates regardless of size 
and temperature feature a BCC Cu wetting layer between the 
BCC Fe matrix and the 9R Cu core. As will be described be- 
low, the latter plays a crucial role in understanding the thermo- 
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FIG. 1. Cross sectional views of Cu precipitates obtained in VC-SGC MC/MD simulations at TOOK exhibiting both (a-c) BCC and (d-f) 
9R structures. The precipitates contain (a) 7,000, (b) 14,900, (c,d) 22,800, and (e,f) 34,600 atoms. The coloring indicates different local 
environments: blue=BCC, red=FCC, green=HCP, white=other. The 9R structure appears as a sequence of two HCP atoms and one FCC atom 
along the [001]9r direction [also compare (g)]. The twin planes are parallel to {llOjg^.^. planes and as can be seen in (f) extend through the 
entire precipitate. Note that the 9R cores of the precipitates shown in (d-f) are enclosed by a BCC Cu wetting layer. The out-of-plane vectors 
in (a-e) are [111]bcc and [IIOJqr, respectively, (g) Schematic of 9R structure adapted from Ref. 4 Triangles, squares, and circles mark atomic 
positions in alternating close-packed planes, top views of which are shown on the left-hand side. The stacking sequence is shown on the 
right-hand side, which demonstrates that if only next nearest neighbors are considered (as in the Ackland- Jones analysis) the 9R unit cell is 
equivalent to three repetitions of the sequence HCP — FCC — HCP. In general, the 9R structure has monoclinic symmetry. For simplicity, here 
we show only the orthorhombic case. 



dynamics of the BCC-9R transition as well as the dislocation- 
precipitate interaction. 

The simulation results are substantiated by the observa- 
tion that the structures of 9R precipitates shown in Fig. [T] 
closely resemble TEM micrographs of 9R Cu-precipitates in 
a-Fe.^ ^ ^ A more detailed analysis of the structure and shape 



of both BCC and 9R precipitates is deferred to Sects. Ill B and 

Next we consider the size and temperature dependence of 
the BCC-9R transition. Figure [2] shows the number of atoms 
in different local environments as a function of precipitate size 
at 700 K. Precipitates with more than about 28,000 atoms un- 
dergo a sharp structural transition as they transform from BCC 
to multiply-twinned 9R structures. 

The simulations also allow us to study the reverse transi- 
tion. To explore this aspect, a series of simulations with de- 
creasing Cu content was carried out starting from a 9R pre- 
cipitate obtained in one of the earlier simulations. The results 
of this investigation are shown by the open symbols in Fig.[2j 
which reveals that the reverse transition occurs for a precipi- 
tate size of approximately 18,000 atoms. There is thus a pro- 
nounced hysteresis between forward (BCC^9R) and back- 
ward (9R^BCC) transition extending from about 18,000 to 
28,000 atoms. This is further illustrated by the BCC and 9R 
precipitates in Figs. [TJc) and (d), respectively, both of which 
contain about 22,800 atoms. It should be stressed that the 
transition can be observed on the time scale of our simula- 
tions because it is martensitic (non-diffusive) in nature and 
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FIG. 2. Fraction of atoms in Cu precipitate in different atomic neigh- 
borhoods at 700 K At smaller sizes BCC motifs dominate indicating 
that Cu precipitates assume the BCC structure. For larger precipitates 
the number of atoms in close-packed environments (9R=FCC-i-HCP) 
is larger as Cu precipitates are now adopting a multiply-twinned 9R 
structure. 



transformations are observed when one structure becomes me- 
chanically unstable with respect to the other. 

This type of analysis was carried out in 100 K intervals be- 
tween 200 and 700 K, the results of which are shown by the 
"horizontal" data points in Fig. [S] Further VCSGC-MC/MD 
simulations were then conducted at constant concentration for 
both decreasing and increasing temperature starting from pre- 
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FIG. 4. Eigenvalues of gyration tensor for Cu precipitates at 700 K. 
The eigenvalues are denoted by Af with A? > A^ > A^. 



FIG. 3. Phase diagram indicating the structure of Cu precipitates 
as a function of temperature and size. The light colored regions 
next to the dark gray phase transition line indicate the extent of the 
hysteresis, which vanishes below 300 K. "Vertical" and "horizontal" 
data points have been obtained from VCSGC-MC/MD simulations 
at constant temperature and concentration, respectively (see Fig.[2|. 
The "MD" data points were obtained from conventional MD sim- 
ulations (without atom type swaps) during cooling and heating, re- 
spectively (see Sect.|niE]). The open diamonds indicate experimental 
data points from Ref.^ 



viously equilibrated BCC and 9R precipitates, respectively. 
The results obtained in this fashion are marked as "vertical" 
data points in Fig. [3] and are found to be consistent with 
the results from sampling at constant temperature. The two 
most important features of Fig. [3] are: (/) the precipitate size 
at which the BCC-9R transitions occurs decreases signifi- 
cantly with temperature and (//) the width of the hysteresis 
between forward and backward transitions systematically de- 
creases upon cooling and vanishes between 200 and 300 K. 

The pronounced hysteresis at high temperatures evident in 
Fig. [3] clearly indicates a first-order phase transition. The van- 
ishing hysteresis with decreasing temperature, however, sug- 
gests that the first-order character diminishes and eventually 
the order of the transition changes. The thermodynamics of 
the transition as well as the strong temperature dependence of 
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FIG. 5. BCC Cu precipitates containing about 2,200 atoms that 
are representative for temperatures of (a) 700 K, (b) 500 K, and (c) 
300 K illustrating the increasing faceting with decreasing tempera- 
ture. The top row shows the atomic positions where blue atoms have 
BCC coordination while red atoms indicate under-coordinated atoms 
at surfaces ledges and corners. The bottom row shows the convex 
hulls. 



where N is the number of atoms in the precipitate, r\ ' denotes 
the component of the position vector of atom i along Cartesian 
direction a and the origin of the coordinate system has been 



the critical transition size are investigated in detail in Sect. IV chosen such that 2_ji 



N J^) 



0. Approximately speaking, the 



B. Structure and shape of BCC precipitates 

As indicated in the previous section, BCC precipitates as- 
sume a spherical shape, while 9R precipitates tend to be ellip- 
soidal. To make this statement more quantitative, the eigen- 
values of the gyration tensor have been computed, which is 
defined as 



1 ^ 



(1) 



gyration tensor describes the ellipsoid that most closely re- 
sembles the point cloud {r^ }. If all eigenvalues are identical 
(A^ = A2 = A3) the point cloud (or in the present case the 
precipitate) is best described by a sphere. On the other hand, 
for a prolate ellipsoid one would obtain one eigenvalue to be 
significantly larger than the other two (A^ > X^ = A3). 

Results of the gyration tensor analysis are shown in Fig. |4] 
for a temperature of 700 K. For BCC precipitates the three 
eigenvalues are identical within the error bars, which is con- 
sistent with the spherical shape observable in Fig.[TJa-c). 

A visual inspection of the precipitates reveals a strong ten- 
dency for faceting as temperature is lowered. While at 700 K 
precipitates are almost perfectly spherical, at 300 K they ex- 



hibit strong faceting toward {110} surfaces with small patches 
of {100} facets 

In Fig.[TJa) time- averaged positions were employed to carry 
out the Ackland- Jones analysis, which reveals that the en- 
tire precipitate has adopted the BCC structure. If instead 
of using time-averaged positions, one employs instantaneous 
positions, FCC and HCP motifs can, however, be found as 
well. This observation is not unexpected since it is well 
known that the phonon dispersion of BCC Cu exhibits un- 
stable modes. ^^^^ MD simulations confirm that bulk BCC Cu 
remains unstable also at elevated temperature and over a wide 
range of volumes. This suggests that the stabilization of BCC 
Cu for precipitates with several thousand atoms is intimately 
connected to their finite size as well as the Fe-Cu interface. 
Indeed, one observes that for thin layers of BCC Cu wedged 
in between BCC Fe the interface tension is sufficient to sta- 
bilize BCC Cu at zero K. However, as the thickness of the 
Cu slab increases, the BCC structure again becomes unsta- 
ble. This implies that the stabilization of BCC Cu precipitates 
with up to 28,000 atoms at 700 K is the result of the combined 
effect of vibrations, finite size and interfaces that individually 
would not be sufficient to stabilize BCC Cu beyond very small 
sizes. 



C. Structure and shape of 9R precipitates 

Next, we analyze the structure of 9R precipitates in the con- 
text of experimental observations and compare the 9R pre- 
cipitates shown in Fig. [TJd-f) with TEM images from e.g., 
Refs. m |5] and [9| In all cases precipitates exhibit the her- 
ringbone pattern that is characteristic of multiply-twinned 
9R structures.^ ^ Both simulated and experimentally observed 
precipitates are notably elongated in the direction parallel to 
the trace of the twin plane with aspect ratios of 1.2 (simu- 
lation) and 1.1-1.3 (experiment). ^^ The orientation relations 
between BCC matrix and 9R domains observed in our sim- 
ulations agree with experiment as well as theoretical predic- 
tions by Kajiwara^^ with the [llOjgR direction being parallel 
to [111] BCC and the (114) qr plane being parallel to (OII)bcc- 
Twinning occurs with respect to {114}9r planes which are 
parallel to {110}bcc planes. 

Based on extensive TEM work, Othen et ah conducted a 
careful analysis of 9R variant thickness (twin spacing). Up 
to a diameter of about 13 nm {i.e., well above the size range 
covered by our simulations) they found twin domains to be 
evenly spaced with a width of about 1.5 ± 0.5 nm. With the 
exception of very small precipitates that contained only two 
segments, they also found a slight increase of twin spacing 
with precipitate size. 

In general, the present simulations agree well with these 
observations. Specifically, we obtain twin spacings between 
0.9 ± 0.1 nm at 300 K and 1.2 ± 0.1 nm at 700 K. Twins are 
evenly spaced and usually contain several structural defects 
such as stacking faults and steps inside the twin planes, very 
similar to the precipitates observed experimentally.^ 

An important structural parameter is the misorientation an- 
gle between {001}9r and {110}bcc planes [see Fig.|6jc) for 



an illustration]. Using the angular deviation of (003) 9r and 
(006) 9R spots from (110)Bcc,Fe spots in their dark field im- 
ages, Othen and co-workers obtained an experimental value 
of 4 ± 0.5° while from fringes in their bright field images they 
extracted a value of 4.3 ± 0.7°. Misorientation angles were 
extracted from the present simulations using time averaged 
atomic coordinates. In the largest precipitates misorientation 
angles in the center of the precipitates are just under 4°, in 
reasonable agreement with experiment. A more detailed look 
at the results, however, reveals a more complex behavior than 
is suggested by the experimental analysis. Figure |6ja) demon- 
strates that misorientation angles vary across the Cu precipi- 
tate. They reach maximal values near the precipitate center 
and drop toward the precipitate-matrix interface. These fea- 
tures are shown in exaggerated fashion in Fig.[6jc). 

The maximum misorientation angle is practically temper- 
ature independent. It does, however, depend on precipitate 
size. Figure |6jb) shows that the maximum misorientation an- 
gle reaches a limiting value for large precipitates (about 3.6° 
for A^ ^ 10^ atoms) and almost completely smoothly ap- 
proaches zero as the precipitate decreases in size. As in the 



case of the internal energy [compare Fig. [TO[a) below] there 
is a small discontinuity left at this temperature. 

We conclude the structural description of 9R precipitates 
with an analysis of their shape in terms of the principal values 
of the gyration tensor as shown in Fig.|4] Unlike BCC precip- 
itates, 9R precipitates display a very pronounced asymmetry 
closely related to the orientation relations that govern BCC- 
9R interfaces. Most noticeably, 9R precipitates exhibit a pro- 
nounced elongation along [110]9r, which is perpendicular to 
the figure plane in Fig. [TJd,e). The other two principal axes 
are [221] 9r and [114]9r with the precipitates being slightly 
longer along the former. As precipitate size increases the as- 
pect ratio between the primary axis and the two shorter axes 
increases. Finally, similar to BCC precipitates, one observes 
a pronounced faceting of 9R precipitates as temperature is re- 
duced with the dominant interfaces being parallel to {110}bcc 
planes. 



D. Transition mechanism from BCC to 9R 

The results presented in the previous sections already pro- 
vide ample evidence for the martensitic character of the BCC- 
9R transition. To complete this analysis, we now describe 
the atomic displacement pattern associated with the BCC-9R 
transition. Figure [7] shows a transformed region in a multiply- 
twinned 9R precipitate projected onto (III)bcc- The atomic 
positions correspond to ideal BCC sites with arrows indicat- 
ing displacements that lead to 9R. The plot reveals a curl-like 
displacement pattern with a periodicity of three atomic layers 
along [001] 9R. The right hand side of Fig. [t] shows the smallest 
unit cell that contains all distinct displacements in a single 9R 
domain. The cell contains 1 8 atoms and is three atomic layers 
thick along the [110]9r direction (perpendicular to the paper 
plane). The atoms also undergo regular displacements along 
[110]9R, which are not visualized here. This path is specific 
for a twin spacing of about 1.2 nm. Other twin spacings arise 
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FIG. 6. (a) Misorientation angle within a single 9R domain parallel to the twin plan for precipitates of various sizes at 300 K. (b) Corresponding 
variation of the maximum misorientation angle with precipitate size. Also shown is the fraction of close-packed atoms in the precipitate. The 
vertical gray bars indicate the precipitate sizes at which dislocation simulations were carried out. (c) Schematic depiction of spatial variation 
of misorientation angle throughout a multiply-twinned 9R precipitate. Misorientation angles are strongly exaggerated for illustration. 
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FIG. 7. Atomic displacement pattern connecting BCC to multiply- 
twinned 9R structure. On the left hand side, a close up view of a 
slab (three atomic layers thick) is shown that is parallel to (III)bcc 
and (110)9R. Filled circles represent ideal BCC positions with the 
arrows indicating displacements that lead to 9R (scaled by a factor 
of three for clarity). Vertical gray bars mark twin planes while the 
slanted gray stripes are parallel to {OOIJqr planes. On the right hand 
side, the smallest unit cell is shown that comprises all distinct dis- 
placements that take place during the BCC-9R transition leading to 
a single 9R domain. Color coding as in Fig.[T] 
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ulation during which the precipitate size was gradually reduced. The 
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multiply-twinned 9R precipitate close to the transition (3,200 atoms), 
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from qualitatively similar curl-like displacement patterns. 

A visual analysis of the simulation results provides the fol- 
lowing picture for the formation of multiply-twinned 9R par- 
ticles. At temperatures at which the transition exhibits clear 
first-order character (^ 400 K), just prior to the transition sin- 
gle 9R domains appear that can extend over significant parts of 
the precipitate along [112]bcc and [111]bcc (compare Fig.[T]) 
but include only a couple of atomic layers along [1 10] bcc- 
These domains are, however, unstable and decay rapidly. In 
contrast, right at the transition, the nucleation of the first 9R 
domain is immediately followed by the nucleation of a second 
9R domain of opposite orientation. This event is crucial for 
stabilizing the first singly-twinned 9R nucleus. Subsequently, 
9R segments are added rapidly as the transformation sweeps 
through the entire precipitate. 

For temperatures ^ 300 K the precipitate size at the tran- 
sition can no longer support two twinned 9R domains. Un- 



der these conditions the transition occurs without hysteresis 
and with vanishing nucleation barrier as shown in Fig. [8] For 
sizes below 2,700 atoms the entire precipitate adopts a BCC 
structure. Between 2,700 and approximately 3,000 atoms the 
precipitate structure fluctuates rapidly between BCC and lo- 
calized close-packed (CP) motifs while the average numbers 
of BCC and CP atoms vary approximately continuously with 
precipitate size. Eventually for larger precipitates multiply- 
twinned 9R structures are stabilized similar to the behavior 
at higher temperatures. Regardless of temperature and size 
9R precipitates exhibit a core- shell structure comprising a CP 
core and a BCC wetting layer, the size of which varies be- 
tween 1.5 and 3.5 atomic layers depending on temperature and 
size. 
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FIG. 9. Fractions of atoms in Cu precipitate in different local en- 
vironments during the course of a MD simulation. During the first 
4 ns of the simulation the temperature decreased linearly from 700 
to 500 K while during the following 4 ns, the temperature was raised 
to its original value. The precipitate contained approximately 18,000 
atoms. The transition temperatures obtained form the present simu- 
lation are plotted as open symbols in Fig. [3] 



E. Quenching simulations 



In Sect. |III A| it was shown that sampling the phase dia- 
gram in Fig. |3] horizontally (constant temperature) and verti- 
cally (constant concentration, approximately constant precipi- 
tate size) yields identical results when using VCSGC-MC/MD 
simulations. Since these simulations involve atom type swaps 
precipitates are allowed to change shape as they pass through 
the BCC-9R transition, which in reality requires diffusion of 
atoms near the Fe-Cu interface. Depending on the cooling 
rate in an experiment diffusion is not always possible and 
shape changes (or their suppression) could affect the phase 
transition. 

To resolve this question a number of conventional MD sim- 
ulations without atom type swaps were carried out, represent- 
ing the limit in which diffusion is entirely suppressed (very 
fast quenching). BCC precipitates of varying size, which had 
been previously equilibrated about lOOK above the expected 
transition temperature, were cooled down by 200 K over a pe- 
riod of 4 ns. Results from a typical simulation run are sum- 
marized in Fig. [9] It is observed that the transition is virtually 
instantaneous, providing further evidence for the martensitic 
character of the transformation. For the precipitate shown in 
Fig. [9] the transformation to a multiply-twinned 9R structure 
occurs at a temperature of about 561 K, which is notably lower 
than the transition temperature of 585 K suggested by Fig.|3] 

To demonstrate that this slight temperature offset is not 
due to the large quenching rate, the temperature was ramped 
back to its original value. The reverse transition occurs again 
slightly below the temperature that one would expect based on 
Fig.|3](633 K v^ 661 K), which implies that the observed tem- 
perature difference cannot be explained by cooling/heating 
rate effects. ^^ 

The temperature shift can rather be rationalized in terms of 
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FIG. 10. (a) Average potential energy of atoms in Cu precipitates as 
a function of precipitate size for different temperatures. The vertical 
black bars represent the heat of transformation, which is plotted in 
(b) as a function of temperature. Lines are guides to the eye. 



the different shapes of BCC and 9R precipitates. In the previ- 
ous sections it was demonstrated that while BCC precipitates 
exhibit a spherical shape 9R precipitates resemble a prolate 
ellipsoid. The MD simulation described above started from 
a spherical BCC precipitate. Since shape changes were sup- 
pressed the spherical shape was also imposed on the 9R struc- 
tures that formed during cooling. The spherical shape being 
less favorable for 9R precipitates increases the meta- stability 
range of BCC precipitates and causes a shift of both forward 
and backward BCC-9R transition to lower temperatures. As 
can be seen in Fig. [3] this behavior is observed consistently 
for precipitate sizes down to about 12,000 atoms. For even 
smaller precipitates the difference in shape between BCC and 
9R precipitates becomes negligible and the shift in the transi- 
tion temperature vanishes. 



IV. THERMODYNAMICS OF BCC-9R TRANSITION 

A. Size and temperature dependence of the transition 

The key result presented hitherto in this paper regards the 
phase diagram in Fig.[3j which exhibits a strong increase in the 
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FIG. 11. Potential energy per atom as a function of distance from precipitate center for BCC and 9R precipitates at (a) TOOK (transition at 
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local environments for a precipitate containing 22,800 atoms at TOOK, (d) Cross section of snapshot of precipitate shown in (c). The wetting 
layer region in (c) appears larger than in (d) due to the implicit projection of the elongated shape onto a sphere. 



size of the BCC-Cu precipitates with temperature before tran- 
sition to the 9R phase occurs. This unusual size-temperature 
phase diagram can be explained by first noting that in bulk, 
the BCC-Cu phase is mechanically unstable. The unstable 
modes of the BCC-Cu phase are stabilized in small precipi- 
tates due to confinement and the restrictions imposed by the 
a-Fe matrix. At finite temperatures, the entropy generated by 
these soft modes, then leads to stabilization of the BCC phase 
relative to the 9R phase. 

Even more interesting is the significant change that is ob- 
served in the character of the phase transformations as a func- 
tion of temperature. At high temperatures, the BCC-9R tran- 
sitions exhibit strong first-order character and prominent hys- 
teresis. Upon cooling the discontinuities become less pro- 
nounced and the transitions occur at ever smaller cluster sizes. 
Below 300 K the hysteresis vanishes while the latent heat di- 
minishes. 

To provide a more quantitative understanding of the char- 
acter of the transition we have analyzed the size and temper- 
ature dependence of the potential energy of the atoms in the 
Cu precipitate^^ as shown in Fig. 10 a). At higher tempera- 
tures the data exhibits a pronounced discontinuity as expected 
for a first-order transition. As the temperature is lowered the 
precipitate size at which the transition is observed decreases, 
and so does the discontinuity in the potential energy. It prac- 
tically vanishes below 300 K, in the same temperature range 
in which the hysteresis disappears. This aspect is illustrated 
in Fig. [TOjb), which shows strong variation of the latent heat 
with temperature. In the remainder of this section we will 
focus on resolving the origin of the unexpectedly large tem- 
perature dependence of the structural transformation of the Cu 
precipitates. 



B. Microscopic analysis 



relation between structure and potential energy profile. BCC 
precipitates exhibit a clearly defined core region, in which the 
potential energy is constant. For a given temperature this en- 
ergy, which can be interpreted as the value equivalent to bulk 
BCC Cu, is size independent.^^ Outside the core radius the 
potential energy drops toward the Fe-Cu interface. This rep- 
resents the wetting layer between BCC Cu core and BCC Fe 
matrix. 

The profiles for larger 9R precipitates [Fig. [TTJa)] also re- 
veal a plateau in the core region, which is, however, smaller 
in terms of radius than for BCC precipitates of the same size. 
Outside this inner core region the potential energy increases 
up to a maximum value ["full core" region in Fig.[TTJc,d)] be- 
fore it decays toward the Fe-Cu interface in a similar fashion 
as for the BCC precipitates ("wetting layer" region). Note that 
the shoulder that appears for example in the wetting layer re- 
gion of the energy profile in Fig.[TTJc) is due to the projection 
of an elongated ellipsoidal precipitate onto a sphere (compare 
Fig. [4]). For the same reason the wetting layer seemingly has 
a larger spatial extent in Fig.[TTJc) than in Fig.[TTJd). 



Figure [TT] shows the potential energy as a function of dis- 
tance from the precipitate center for both BCC and 9R pre- 
cipitates at two different temperatures, and illustrates the cor- 



As shown in Fig.[TT]the different regions in the energy pro- 
file can be correlated with different atomic structures. The 
inner core region corresponds to the close-packed core of the 
precipitate representing bulk 9R, whereas the outermost re- 
gion comprises BCC atoms and acts as a wetting layer be- 
tween Cu and Fe. The intermediate (outer core) region is char- 
acterized by a continuous variation in the average number of 
BCC and close-packed atoms as well as the potential energy. 
This region is required to connect the BCC and close-packed 
lattices crystallographically, and corresponds to the outermost 
layers of the 9R region of the precipitate for which the misori- 
entation angle is reduced from its maximum value as shown in 
Fig.|6ja). It thus corresponds to a close-packed region under 
shear strain. As will be discussed below this region is crucial 
for understanding the strong reduction in the latent heat with 
decreasing temperature (and transition size). 

For large precipitates (^ 7, 000 atoms) the inner core 
("bulk") potential energy of 9R precipitates is practically size 
independent as exemplified in Fig.[TTJa). For smaller precip- 
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itates, however, the inner core region can no longer fully re- 
lax since the confinement constrains the misorientation angle, 
which raises the energy compared to larger 9R precipitates 
[Fig.llTJb)]. 

We now discuss in detail the core and wetting layer contri- 
butions to the latent heat for different temperatures as shown 
Fig. [12] Both contributions approach asymptotic values, 
which for the core amounts to 20meV/atom. This is about 



half the value of 40 meV/atom that is found for the energy 
difference between the ideal BCC and 9R structures at zero 
Kelvin. The lower value obtained from direct analysis of the 
precipitate core energies can be attributed to the high planar 
fault density in actual 9R precipitates [twin planes as well as 
stacking faults, see Figs. [T] and [TTJd)]. 

While the core potential energy difference between BCC 
and 9R precipitates is positive. Fig. [12] shows that the wetting 
layer energy difference has the opposite sign. This feature, 
in addition to the wetting layer increasing in size in smaller 
precipitates, is responsible for the drastic reduction in total 



latent heat as temperature is lowered (see Figs. [TO] and [12]). 
This is because at lower temperatures the transition occurs at 
smaller precipitate sizes, where the number of Cu particles 
belonging to core and wetting layer begin to equalize, leading 
to mutual cancellation of their contributions to the latent heat. 
We conclude the description of energy profiles by provid- 
ing further insight on the effect of strain on the core energet- 
ics, which resides in the misorientation angle data presented 
in Fig. [6] and Sect. IIIC Figure ^2i) demonstrates that the 
misorientation angle is not constant across the precipitate but 
rather increases continuously from zero near the interface to 
its maximum value in the core region of the precipitate. This 
observation explains the gradual change in the potential en- 
ergy profile in the outer core region in Fig. [TT] As already 



mentioned in Sect. IIIC the maximum misorientation angle is 
practically temperature independent and levels off at approxi- 
mately 3.6° in the limit of large precipitates. It is fundamen- 
tally determined by the 9R twin boundary energetics. Using 
the data in Fig.[6];b) one can demonstrate a simple correlation 
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between misorientation angle and the core contribution to the 
latent heat. Figure [13] shows a fit to the data from Fig. ^h) 
and indicates the maximum misorientation angles that are ob- 
served at transition for different temperatures. While the data 
points above 500 K all yield approximately the same maxi- 
mum misorientation angle, the 300 and 400 K data points cor- 
respond to considerably smaller values. As shown in the inset 
the core potential energy difference exhibits a simple linear 
dependence on the maximum misorientation angle and is ac- 
cordingly correlated with strain. 
As can be seen in Fig. 



12 the confinement strain has a 



stronger effect on the energy difference in the core region than 
in the wetting layer. In fact below 300 K when the transi- 
tion occurs without hysteresis, the misorientation angle also 
becomes vanishingly small. According to the fit in Fig. 13 
the latter should become zero at a size of about 2,800 atoms, 
which correlates very well with the size range between 2,700 
and 3,000 atoms, for which one can no longer observe stable 
9R twins in simulations at 200 K (see Fig.[8]and Sect. HID ). 



V. DISLOCATION-PRECIPITATE INTERACTION 

Ultimately, Cu precipitation in ferritic alloys leads to hard- 
ening by impeding dislocation motion. Hardening is a macro- 
scopic quantity that takes into account multiple dislocation- 
precipitate interactions over a statistically significant volume. 
However, the elementary hardening mechanism in our case 
is dislocation bowing around a precipitate, which occurs at 
the atomic scale and is amenable to study by MD simula- 
tions. In fact, recent simulations show that the dislocation 
core, which represents the very atomic essence of a disloca- 
tion, may su ffer st ructural transformations when near or inside 
a precipitate .2233 

Similarly, in the presence of a stress field, such as that of 
a dislocation, grain boundary, etc., Cu precipitates may suffer 
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internal changes. Indeed, Cu BCC precipitates have been ob- 
served to locally transform to a close-packed structure when 
subjected to the stress field of a screw dislocationPJEol 

The incremental increase in shear stress r^ due to Cu pre- 
cipitates can be estimated from a dispersed barrier hardening 
(DBH) modeP^ with knowledge of size, number density and 
composition of precipitates, "^"^ 



~17 



(2) 



where fi is the matrix (BCC Fe) shear modulus, h is the Burg- 
ers vector, L is the average barrier spacing, and a is the 
strengthening coefficient. Estimates of a can be obtained by 
recourse to line tension expressions of the DBH, such as the 
Fleischer-Friedel formula,!^ 
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where ^c is the critical cutting angle after which the disloca- 
tion is able to traverse the obstacle. Measurements of (^c can 
be obtained via, e.g., in-situ TEM observations of dislocation- 
irradiation obstacle interactions."^^ MD can provide direct es- 
timates of Tp, although it suffers from length and time scale 
limitations that result in unphysically high strain rates. ^^E^^EH 

The crystallographic orientation of the precipitate adds yet 
another dimension to the already multidimensional character 
of the dislocation-precipitate interactions. Most of the pub- 
lished works on dislocation-precipitate interactions in Fe-Cu 
alloys consider small coherent BCC precipitates,^^ "^^ ^^ both 
in order to study incipient hardening at the beginning of Cu 
precipitation and due to the absence of high-fidelity atom- 
istic models of 9R and FCC Cu precipitates. Our work ef- 
fectively eliminates this limitation for the first time. Here, 
we study screw dislocation interactions with several precipi- 
tates obtained by quenching from 600 to 10 K including both 
BCC and 9R configurations, and compute Tp as a function of 
L. Although statistically speaking the average barrier spacing 
in Eq. O L = {2Nrc^)~ ' is taken from the known num- 
ber density TV = (84a x 306 x 52c) ~^ and diameter D of the 
precipitates, here we simply take L = 84a — D for consis- 
tency with the MD setup. We then measure the critical shear 
stress required to overcome the obstacle, the results for which 
are given as Tp-L data pairs in Table [l| Snapshots from a typ- 
ical simulation illustrating the process of a dislocation cutting 
the precipitate are shown in Fig. [14] 

Data points corresponding to the [llOjgR orientation can 
be fit to an equation of the form Tp = AL~^, which yields 
A = 8.42 X 10"^ GPanm. Using valueJ^of 6 = 0.25 nm 
and /i = 71.5 GPa, we obtain a = 0.47 for the strengthening 
coefficient in Eq. ([2]), which using Eq. ([3| results in an effec- 
tive exit angle of 105°. This is in good agreement with the 



direct analysis of the simulations as shown in Fig.[T4[d). 

Regarding the interaction mechanism, in the previous sec- 
tions it has been demonstrated that the (111) directions of the 
Fe matrix are aligned with the (110) directions of 9R pre- 
cipitates, which are close-packed. This implies that when 
a ^(111) screw dislocation line is co-linear with this close- 
packed direction, the screw character of the dislocation can be 



TABLE L Calculated precipitate strength for precipitates of differ- 
ent sizes where L is the dislocation line length and D the precipitate 
diameter parallel to the dislocation line. The orientation column indi- 
cates the vector in the frame of reference of the precipitate that aligns 
with the screw direction in the matrix. The last precipitates are nom- 
inally equal in size but the orientation with respect to the dislocation 
line differs. As a result of the elongated shape of the 9R precipitate 
(compare Figs. [T] and |4]) the latter has a shorter effective diameter 
projected onto the dislocation line. 
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conserved within the precipitate. Then, the resistance exerted 
by the precipitate on the dislocation originates from the energy 
cost associated with the core energy of the 9R dislocation as 
well as from the herring-bone slip plane structure encountered 
inside the precipitate. This structure must be negotiated via 
multiple cross-slip steps that may hinder dislocation motion 
as well. 

When the principal (110)9r axis is not aligned with the 
screw direction, the dislocation must nucleate a set of inter- 
face dislocations to conserve the Burgers vector. In that case, 
the energetics associated with these nucleation episodes result 
in considerably higher hardening. To exemplify this effect, 
we have studied one such interaction corresponding to the ori- 
entation [114]9R, which results in a critical stress that is two 
times higher than for the first orientation. The calculated Tp 
value is also given in Tab. [T| and included for comparison in 
Fig.fBl 



VI. DISCUSSION 

A. Thermodynamic implications 

The transition diagram in Fig. [3] represents the first key re- 
sult of the present study. Most importantly it demonstrates 
that the precipitate size at which the transition occurs varies 
strongly with temperature, changing from 2,700 atoms (equiv- 
alent to a precipitate diameter of '^ 4 nm) at 200 K to 22,800 
atoms (~ 8nm) at 700 K. This strong size dependence is a 
direct result of the soft vibrational modes of the BCC Cu 
phase, which are responsible for this phase having a larger 
vibrational entropy than 9R Cu. Hence increasing tempera- 
ture leads to stabilization of larger BCC-Cu precipitates. As 



will be discussed in Sect. VI B below the temperature depen- 
dence of the transition size is in very good agreement with the 
interpretation of experimental data obtained by Monzen and 
co-workers^. 

The second important finding pertains to the order of the 
transition. At higher temperatures it exhibits pronounced 
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FIG. 14. Illustration of screw dislocation (line vector parallel to [111]bcc and [IIOJqr) interacting with a 9R precipitates containing 17,400 
atoms (D = 8.8 nm, compare Tab.|l]). The screw dislocation moves along [110]bcc as indicated by the dashed line, (a) Initial precipitate, (b-d) 
Dislocation inside precipitate just prior to detachment, (e-h) Dislocation as it exits precipitate. The magnitude of the exit angle expected based 
on Eq. ^ is indicated in (d). 



hysteresis and a finite latent heat (6 meV/atom at 700 K, see 



Fig. 12). With decreasing temperature the latent heat dimin- 
ishes and the nucleation barrier vanishes at approximately 
300 K, whereupon the transition occurs without hysteresis. 
The disappearance of the nucleation barrier is accompanied 
by characteristic changes in the misorientation angle. For pre- 
cipitates below approximately 6,000 atoms (5 nm) it quickly 
decreases and reaches zero at a size of 2,800 atoms, which 
correlates with the critical size observed at 200 K. The strong 
size dependence of the misorientation angle is the result of 
the crystallographic boundary conditions along the interface 
between close-packed (inner core in Fig. [TT]) and BCC Cu 
(wetting layer) regions, which are manifested in the presence 
of a layer of close-packed Cu atoms under shear strain (outer 
core). 

The ingredients that underlie the strong size dependence 
and the change in transition character are: (/) two crystallo- 
graphically mismatched boundary phases and (//) the minor- 
ity component is dynamically unstable or at least vibrationally 
soft in the lattice structure of the matrix phase. These condi- 
tions are very general and can be found in a large number of 
systems. The features described in the present work should 
therefore also be observable in other systems. 

Vacancies are known to pla y an important role in the forma- 
tion of Cu precipitates in iron-lS^SEll They have also been dis- 
cussed in the context of the BCC-9R transformation (see e.g., 
Ref.[T3land[T6l). While thermodynamically vacancies are not 
needed for the transition they are likely to play a crucial role 
for the kinetics of the transformation and the loss of coherency 
between BCC matrix and 9R precipitates (see Sect. | VI B] ). 



B. Implications for interpretation of experiments 

Most TEM studies (see e.g., Refs.|4l|5l and[9l) were carried 
out at room temperature on samples that were quenched from 
temperatures above 750 K to room temperature. The smallest 



9R precipitates observed in these studies were typically just 
slightly larger than 4 nm in excellent agreement with the crit- 
ical size at about 300 K obtained in the present simulations. 

In the present context comparison with the experimental 
study by Monzen and co-workers^ is particularly interesting. 
Motivated by some earlier results by Habibi-Bajguirani and 
Jenkins^, these authors carefully investigated the size depen- 
dence of the angle Lp enclosed by {OOIJqr planes in adjacent 
twin segments in 9R precipitates as a function of annealing 
temperature. 

They observed that precipitates annealed for one hour at 
673 K exhibit a value of (^ ^ 128° up to a size of about 9 nm. 
For larger precipitates the angle drops sharply to (/p ^ 122°. 
The size at which this transition takes place is strongly depen- 
dent on the annealing temperature, with lower temperatures 
leading to the transition occurring at smaller sizes. Monzen et 
al interpreted this finding as an indication for a temperature- 
dependent transition between BCC and 9R structures strongly 
reminiscent of the present findings. The experimental data 
points are included in Fig. [3] and in fact display a remarkable 
agreement with the BCC-9R transition line obtained from 
simulation.To appreciate the interpretation by Monzen et al 
and the connection to the present simulations let us revisit 
the experimental observations in more detail. As discussed 
extensively in Sect. |IIIC| {QQIJqr planes are approximately 
parallel to IIObcc planes. The deviation from perfect align- 
ment is measured by the misorientation angle (see Fig.|6]). The 
value of (/9 = 128° observed by Monzen and co-workers for 
smaller precipitates is equivalent to a misorientation angle of 
about ^ = 4° , which agrees both with earlier measurements 
on quenched precipitates'^ ^ ^ and our simulations. The value 
oiLp = 122° obtained for larger precipitates on the other hand 
corresponds to a misorientation angle i3 of only 1 ° . The latter 
is the result of — presumably diffusion mediated — reshuffling 
at the interface, by which the strain in the 9R precipitates is 
reduced.^ It is likely that vacancies play a crucial role in this 
process (compare Refs..l3.and.l6j. 
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A precipitate that adopts the 9R structure at the anneaUng 
temperature will be able to rearrange its interface during the 
extended annealing period and thereby reduce its misorienta- 
tion angle. In contrast a precipitate that adopts the BCC struc- 
ture under the annealing conditions will rapidly pass through 
the BCC-9R transition during quenching with insufficient 
time and activation to transform its interface. As a result it will 
exhibit a misorientation angle of about 4° as dictated by the 
interface alignment with the Fe BCC matrix and the Cu BCC 
wetting layer. This distinction enables the indirect experimen- 
tal detection of the BCC-9R transition by TEM observation at 
room temperature. 

In our simulations we directly observe the diffusion-less 
transition that takes place during rapid quenching (or heat- 
ing). The time scale of our simulations is too short to allow 
for a diffusional transformation of the interface and as a re- 
sult the simulated precipitates exhibit maximum misorienta- 
tion angles of up to 4° . It is important to stress that the experi- 
mental observation is only possible because of diffusion along 
the interface since otherwise 9R precipitates that undergo a 
transition during quenching would be indistinguishable from 
precipitates that adopt the 9R structure already during anneal- 
ing. In this sense experiment and simulation provide compli- 
mentary views, and while the direct simulation of the BCC- 
9R transition confounds the interpretation of the experimental 
observation, the good agreement with experiment verifies the 
simulation. 

It remains to point out that a diffusional interface transfor- 
mation is unlikely to be of importance for the low temperature 



1.5 



behavior described in Sect. |IV] not only because diffusion is 
hampered at such low temperatures but — more importantly — 
because for the relevant precipitate sizes the misorientation 
angle is already strongly reduced due to confinement (com- 
pare Fig. |6]). 

To conclude this section a few additional comments regard- 
ing the structure and shape of Cu precipitates are in order. In 
Sects. IIIB and IIIC it was found that BCC precipitates adopt 
a spherical equilibrium shape whereas 9R precipitates resem- 
ble an elongated ellipsoid. Conversely, in TEM investigations 
smaller 9R precipitates appear often to be characterized as cir- 
cular (in 2-D projection), whereas larger 9R precipitates show 
an elongated shape. This observation is readily understood 
on the basis of the present data as smaller circular 9R pre- 
cipitates seen in TEM investigations actually correspond to 
(spherical) BCC precipitates at annealing temperatures. Dur- 
ing rapid quenching to room temperature these precipitates 
transform into the 9R structure. Shape change, however, is 
suppressed since it once again would require diffusion. In 
contrast larger precipitates adopt the 9R structure also under 
annealing conditions and remain 9R during cooling. This sit- 
uation resembles the MD simulation shown in Fig. [9] 

The simulated precipitates exhibit some features that are ei- 
ther impossible or very difficult to detect experimentally. For 
example Fig.[TJb) reveals that multiply-twinned 9R Cu precip- 
itates have a BCC Cu shell that typically comprises about two 
to three atomic layers and serves as a wetting layer between 
the close-packed Cu core and the BCC Fe matrix. There is 
thus no direct 9R Cu and BCC Fe interface, which is of ut- 
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FIG. 15. Critical shear stress as a function of harmonic average of 
L and D for MD simulations of screw dislocation-precipitate inter- 
actions performed here. Results are normalized to /nb/L to remove 
the dependence with L and /x and facilitate comparison with earlier 
MD simulations."^ 



termost importance for the thermodynamics of the transition 
as discussed in detail above. The existence of a wetting layer 
has some additional implications including the following: (/) 
Experimentally, the diameter of Cu precipitates is typically 
measured using contrast produced by different lattice struc- 
tures. As this probe is insensitive to the BCC Cu shell it leads 
to a slight underestimation of precipitate size. (//) A disloca- 
tion interacting with a 9R Cu precipitate will not see a sharp 
BCC Fe — 9R Cu interface but instead pass pass through three 
distinct zones (BCC Fe, BCC Cu, and 9R Cu). 



C. Dislocations 

In Sect. [V] it was shown that Cu precipitates can be a sub- 
stantial source of resistance to screw dislocation motion in 
Fe-Cu alloys, and that for a given precipitate structure and 
orientation, hardening follows a DBH-type law. For that par- 
ticular structure/orientation combination, the screw direction 
is a principal axis of the precipitate and there is no kinematic 
constraint to the dislocation traversing the obstacle; only the 
energetics of generating a dislocation segment in a different 
crystal structure are responsible for the observed hardening. 
Since we have obtained the hardening response from five pre- 
cipitate sizes, including three different orientations, it is more 
informative to represent the values given in Tab. |l] as a func- 
tion of the harmonic mean between L and D. This is done in 



Fig. 15 which shows normalized critical stresses. 

The figure reveals structure and orientation dependencies 
by placing different structures and/or orientations within dis- 
tinct classes of behavior, as the dashed lines intended as guides 
to the eye reveal. Moreover, as the figure illustrates, our 
BCC Cu point is in excellent agreement with previous cal- 
culations by using a different interatomic potential Bacon and 
Osetsky,^^ providing indirect verification for the dislocation- 
9R precipitate calculations. 
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As another useful verification exercise, we can confirm the 
interaction mechanism based on dislocation renucleation in- 
side the precipitate (enabled by orientations alignment). To 
this end, more sophisticated models, such as the Russell- 
Brown modulus interaction model,^^ that include the modu- 
lus mismatch between the matrix and the precipitate, must be 
invoked. 



r„ = 






E 



(4) 



where /i, 6, and L are defined as in Sect.jVJ and Eq^ and E are 
energies per unit length of the dislocation in precipitate and 
matrix, respectively. For screw dislocations, the ratio Eq^/E 
is equivalent to^^ 



E ~ log(roo/ro) 



log 



^0 



log 



^Cu 



(5) 



where r^ and Too are the dislocation cutoff radius and outer 
cutoff radius used to compute the dislocation energy, /icu is 
the shear modulus of the precipitate (which we approximate 
as that of FCC-Cu), and rcu is the precipitate radius. 

Equation ^ can be fitted to data from the interatomic po- 
tentials used in this work: ji = 71.5 GPa, jiqu = 41.6 GPa 
(again, rotated to the slip geometry used here from the nomi- 
nal value of 76 GPa), r^ = 26, r^o = lOOOro, and rcu = 46. 
This yields a value of a = (1 — Equ/E) ' ^ 0.5, in good 
agreement with the MD analysis. This confirms that disloca- 
tion nucleation inside the precipitate is the dominant source of 
hardening, which is not surprising in light of the fact that we 
have chosen orientation relations that favor this mechanism 
over nucleation of interface partials. In any case, this obser- 
vation runs counter to the conclusion extracted by Fine and 
Isheim^^ in their review of hardening models. Note that the 
values of ro, Too, and rcu are somewhat arbitrary, although 
due to their logarithmic dependence their effect is likely to be 
small. We note that our MD values are considerably smaller 
than those obtained using dislocation dynamics simulations, 
although these contain only a crude model of precipitate struc- 
ture that may not yield the correct hardening .1^3 

In this context the observation of dislocation induced 
martensitic transformations in BCC Cu precipitates^^ ^^ that 
was alluded to above is also of interes t.!^ In MD simulations 
carried out at lOK Shi and co-workers^^ESl observed the par- 
tial and largely reversible transformation of Cu atoms into a 
close-packed structure 9R-like structure in the presence of a 
screw dislocation for BCC precipitates with a diameter of at 
least 2.5 nm (about 690 atoms). The fraction of transformed 
atoms was found to be strongly dependent on the distance 
between dislocation and precipitate surface. For distances 
above 1 .5 nm no transformation was observed up to the largest 
precipitate size considered (4nm equivalent to about 2,800 
atoms), while for small distances more than 30% of the Cu 
atoms transformed to a close-packed structure. 

Bacon and Osetsk>^^ carried out a similar analysis as a 
function of temperature albeit for larger precipitates (6 nm, 
about 9,700 atoms). In their simulations more than 50% of the 
Cu atoms underwent a transformation at about 100 K while 



this fraction dropped to about 10% at 600 K. These data are 
consistent with the results reported in the present paper, which 
provide additional insight into the origin of the observed be- 
havior. The analysis of the temperature dependence of the 
BCC-9R transition (Sect. [Ill]) and its energetics (Sect. |IV]) 
provides a rationale for the strong temperature dependence 
observed by Bacon and Osetsky. It also resonates with their 
conclusion that the temperature dependence of the interaction 
of dislocations with Cu precipitates in a-iron is more complex 
than for voids in the same material. ^^ 

The correlation between the stress field of dislocations and 
the martensitic transformation in BCC Cu precipitates as well 
as the well-known shear instability of BCC Cu^^ ^^ suggest 
that the phase line between BCC and 9R Cu might not only be 
strongly temperature but also shear stress dependent. 



VII. CONCLUSIONS 

We have studied Cu precipitates in Fe-Cu alloys using 
atomic scale simulations and empirical potentials. In terms 
of quantitative information, the main findings of this paper 
are (/) the polymorphic phase diagram of Cu precipitates in 
BCC Fe, and (//) the strength (resistance to dislocation mo- 
tion) of Cu precipitates with different structure. Finding (/) is 
validated by experimental data found in the literature (which 
provide shape, structure, and thermal information), whereas 
(//) is verified against published MD simulations. 

The BCC-9R transition is found to be strongly size and 
temperature dependent. This is explained in terms of the en- 
tropic stabilization of BCC versus close-packed Cu. While 
at high temperatures the transition is first-order it takes 
place without hysteresis at temperatures below approximately 
300 K. The latter behavior correlates with the maximum mis- 
orientation angle, which can be used to estimate the size of 
the nucleation barrier for the structural transformation. 

In addition we have extensively characterized the structure, 
morphology and strengthening effect of Cu precipitates. The 
present approach could conceivably be employed to deter- 
mine the alloy composition and thermal treatments to tune the 
amount of precipitation hardening desired in an alloy. 
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